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The quantum Rabi model describes the ultrastrong interaction of a two-level atom coupled to a 
single quantized bosonic mode. As compared to the Jaynes-Cummings model, in the Rabi model the 
absorption and emission processes do not need to satisfy energy conservation and the usual rotating 
wave approximation (RWA) breaks down. As a result, the atom-field dynamics in the Hilbert space 
splits into two independent parity chains, exhibiting a collapse-revival pattern and exact periodic 
dynamics in the limit of degenerate atomic levels. Here we introduce a mixed Rabi Jaynes-Cummings 
model by considering a three-level atom interacting with two quantized bosonic fields, in which the 
RWA is made for one transition (with a weak atom-field coupling) but not for the other one (with 
an ultrastrong atom-field coupling). As a result, we show that the field in the weak coupled atomic 
transition can be used as a tool to control the atom-field dynamics of the other (strong coupled) 
transition, thus realizing an effective two-level quantum Rabi model with a controllable field. In 
particular, a periodic temporal dynamics of the atom-field state can be realized by appropriate 
tuning of the weak control field, even for non-degenerate atomic levels. A photonic simulator of 
the mixed Rabi Jaynes-Cummings model, based on light transport in evanescently-coupled optical 
waveguide lattices, is also briefly discussed. 


I. INTRODUCTION 

The well-known quantum Rabi model [H-Q , describing 
a two-level atom coupled to a quantum harmonic oscil¬ 
lator, continues to produce rich and surprising physics, 
with plenty of applications in a variety of physical sys¬ 
tems. The quantum Rabi model has been applied to nu¬ 
merous experimental systems in quantum optics or con¬ 
densed matter, such as cavity quantum electrodynamics 
^ED) quantum dots [p , superconducting qubits 

and trapped ions (ill. Hp. In most cases, when 
the external field is weak enough [I^, the rotating wave 
approximation (RWA) is applied and in such way the 
famous Jaynes-Cummings model is obtained [l^. How- 
ever, in recent years, new regimes have been explored 
(iSl - l^ , in which the effect of counter-rotating terms can¬ 
not be neglected. Such regimes are the ultrastrong cou¬ 
pling of light-matter interactions [iTl - f^ and the deep 
strong coupling (DSC) In the DSC regime, the 

absorption and emission processes do not need to sat¬ 
isfy energy conservation and the atom-field dynamics is 
more involved and splits into two parity chains in Hilbert 
space. As a result, the atom-field state undergoes revival 
and collapse dynamics in Hilbert space . Remarkably, 
in the limit of degenerate atomic levels the dynamics be¬ 
comes exactly periodic [l^. Recent works have shown 
that the quantum Rabi model can be simulated by us- 
ing light transport in engineered waveguide superlattices 
2J-[^ . This could allow the DSC regime, which is hard 
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to access experimentally in cavity QED, to be success¬ 
fully simulated in other physical contexts. In spite of 
the vast research in this area and the relative simplicity 
of the Rabi model, its integrability has been proven just 
recently [T7I - |2^ . 

In the past few decades, several theoretical and experi¬ 
mental works have shown that many interesting coherent 
phenomena can be observed when more than two atomic 
levels are involved in the dynamics. In particular, the 
three level system exhibits a plethora of coherent phe¬ 
nomena such as two-photon coherence [s^ , resonance Ra¬ 
man scattering [sij , double resonance proces s , three- 
level super radiance and quantum jumps |34||. to men¬ 
tion a few. The quantum dynamics of a three-level atom 
interacting with two resonant or near resonant modes of a 
steady field has been generally studied by assuming either 
quantized fields, using a dressed-state formalism within 
the RWA (extended Jaynes-Cummings models [^-[37|). 
or the Floquet approach for classical fields [38l - l4j]| . In 
such previous works, the strong coupling regime was con¬ 
sidered solely in the semiclassical limit, i.e. for classical 
fields, where a Floquet analysis of the underlying time- 
periodic equations for the atomic population amplitudes 
can be employed for sinusoidal external fields . How¬ 
ever, for quantized fields the ultrastrong coupling regime 
in the three-level atomic system, which breaks the RWA, 
was not investigated. 

In this paper, we study in details a three-level atomic 
system interacting with two quantized fields, where one 
field is near-resonant and weakly coupled with one atomic 
transition whereas the other field is strongly coupled to 
the other atomic transition. Such a quantum model can 
be refereed to as mixed Rabi Jaynes-Cummings model, 
because the RWA can be applied to one transition, but 
not to the other one. We derive the coupled differential 
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FIG. 1 . Schematic of a three-level system interacting with two 
bosonic modes with n quanta at frequency cJi and k quanta 
at frequency uj2. 

equations, describing the temporal evolution of the quan¬ 
tum system in Hilbert space, and show that the weak- 
coupling filed can be used as a tool to control the dy¬ 
namics in Hilbert space of the atom-field state for the 
other transition. In particular, the weak control field can 
be tuned to realize exact periodic of the atom-field state 
even if the strongly-coupled atomic levels are not degen¬ 
erate. A possible physical implementation of the mixed 
Rabi Jaynes-Cummings model, using arrays of coupled 
optical waveguides with engineered coupling constants, 
is also briefly discussed. 

II. THE MODEL 
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FIG. 2 . Goupling scheme after RWA for the | 2 )at — | 3 )at tran¬ 
sition. Depending on the initial condition, the chain contains 
only (top) even number or (bottom) odd number of bosons in 
the first mode when the atom is in state |l}at. 

performed for one of the two transitions (like in standard 
Jaynes-Cummings model) but not for the other one (like 
in the quantum Rabi model). To study the exact tempo¬ 
ral evolution of the atom-field state |'I'(t)) in the mixed 
Rabi Jaynes-Cummings model, let us expand the state 
vector of the system as 


We consider a three-level quantum system, with atomic 
states |l)atj | 2 )atj and | 3 )at, interacting with two bosonic 
(e.g. electromagnetic) helds of states |n, fc)bos, where n 
and k are the number of bosons in the two fields (Fig.[T]). 
Such system is described by the Hamiltonian 

3 

H = EiGii -|- hijJia\ai + fyj02a\a2 (1) 

i=l 

+ %i(cti 3 + o-3i)(ai -I- a|) -I- hg2{a23 + cr32)(a2 -b 4), 


= 

E [<iW|l)at +<lW|2)at +Cji(t)|3)at 


1 ^ 5 ^)bos j 


(3) 


where is the probability amplitude to have (n, k) 
bosons in the two fields and the atom in level |/)at- Sub¬ 
stitution of the ansatz ([S]) into the Schrbdinger equation 

= H\^{t)), ( 4 ) 


where wi_2 are the frequencies of the fields that are re¬ 
sponsible for the |I)at - |3)at and |2)at - |3)at atomic 
transitions, respectively. The coupling strengths are pa¬ 
rameterized by gi and g2, a\ 2 and 01,2 are the cre¬ 
ation and annihilation operators for the two modes and 
CTij = I*)at at(j|- We now assume <C 0 J 2 and a resonant 
and weak coupling for the |2)at — |3)at transition, hence 

= E3 — E2 and hg2Vk <C £’3 — ^2. Under these 
conditions we can apply the RWA for this transition, dis¬ 
carding the counter-rotating terms 173202 and 172302, and 
the Hamiltonian now reads 

3 

H = Eiau + huJia\ai + hjjJ 2 a\a 2 ( 2 ) 

i=l 

+ %i(oi3 + cr3i)(ai + oj) -b %2(o32a2 + 0230 ^). 

The model described by the Hamiltonian ([ 2 ]) is a mixed 
Rabi and Jaynes-Cummings model, because the RWA is 


yields the following coupled differential equations for the 
probability amplitudes, 

+ {nhuj3 + ktu.2) ( 5 a) 

-b -b 

= E 2 C^nl + + fcSc. 2 ) Ug ( 5 b) 

+ % 2 V^C'S_, 


+ {nh^i + kh.2) cS ( 5 c) 
-b -b 

-b % 2 V'fc-blC^Ei ■ 


The coupling between the amplitudes can be visualized as 
two sets of uncoupled chains, where each lattice site cor¬ 
responds to a different state of the atom-field system. De¬ 
pending on the initial condition, one of the two sets is re¬ 
alised and the other one is irrelevant. The chains in each 
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set are uncoupled in the k direction, because of the RWA, 
while in the n direction they are semi-infinite with a gra¬ 
dient in the boson number. This is depicted in Fig. [H 
where we show the two coupling schemes, for a particu¬ 
lar value of the second-mode boson number k = k. The 
top part of the picture shows the coupling scheme with 
only even number n of bosons in the first mode when the 
atom is in state |l)at- On the opposite, the bottom part 
contains only odd number of bosons n when the atom is 
in |l)at- In what follows, we will assume that the initial 
condition is |'I'(0)) = |l)at|0, K)bos and hence the top part 
of Fig. Ois relevant. A more suitable basis for the vec¬ 
tor state of the system is [|6i), I62), I63), 1^4), • •], 


where 

l^l) = |l)at|0,K)bos, 

1 ^ 2 ) = |3)at|l,K)bos, 

1 ^ 3 ) = |2)at|l, K, + l)bos, 

I 64 ) = |l)at|2,K;)bos, ( 6 ) 

l^s) = |3)at|3,/t)bos, 
l^e) = |2)at|3, K + l)bos, 


In this basis, the Hamiltonian can be represented in the 
following matrix form: 


H = 


El hgiVl 0 

hgi^/l E^ + huji hg2\fK 

0 hg 2 ^/K E 3 + huJi 

0 hgl ^/2 0 

0 0 0 

0 0 0 

0 0 0 


0 

hgiy /2 

0 

E\ -\- 2 ,huji 
hgiVi 
0 

0 


0 

0 

0 

51^3 
+ Shuji 
hg2y/K 

hgi'/i 


0 

0 

0 

0 

hg 2 ^/K 
Eh 3 huji 

0 


0 

0 

0 

0 

hgiVl 

0 


(7) 


where we have discarded the common diagonal term Khuj2 
and we have used the resonance condition E^ = E2 + fiuj2 ■ 
It should be noted that the type of Hamiltonian matrix 
given by Eq. obtained by projection of the full quan¬ 
tized Hamiltonian in the atom-field basis |^)at |n., fc)bos 
{I = 1 , 2 , 3 ), has a similar structure of the Floquet ma¬ 
trix that one would obtain in the semiclassical limit, i.e. 
when the strong field at frequency wi is assumed to be a 
classical field (coherent state). In the semiclassical anal¬ 
ysis, the equations of the atomic field amplitudes with¬ 
out the RWA for the transition | 2 )at — | 3 )at are coupled 
by the oscillating field at frequency wi, i.e. they form 
a set of non-autonomous equations with periodic coef¬ 
ficients. The semiclassical time-periodic Hamiltonian of 
the atomic amplitude probabilities can be cast into an 
equivalent infinite-dimensional time-independent matrix 
Hamiltonian by standard Floquet approach [d^. The 
resulting Floquet matrix has a form similar to Eq. 0. 
The main difference between the classical and quantized 
field descriptions is that in the latter case the couplings, 
described by the two second main diagonals of Eq. 
depend on the photon number. Since the spectrum of 
the matrix H is generally not equally spaced, the atom- 
field dynamics is generally aperiodic, regardless the field 
at frequency uji is considered classical or quantized. In 
particular, if the control field at frequency UJ2 is switched 
off {g2 = 0 ), the Hamiltonian El reduces to that of the 
two-level quantum Rabi. In this case it is known that 
the spectrum of H is not equally-spaced, unless the de¬ 
generate level limit E3 — —>• 0 is considered [l^. In 


the semiclassical limit, the lack of periodicity can be ex¬ 
plained on the basis of Floquet theory of the strongly- 
driven two-level atomic transition, in which two gener¬ 
ally incommensurate frequencies are involved (the fre¬ 
quency uii of the external field and the difference of the 
quasi energies; see, for instance, 0). To calculate the 
energy spectrum of El and the atom-field dynamics in 
the general case g2 0 , one has to resort rather gener¬ 
ally to a numerical analysis. However, some interesting 
physics can be found if one takes a different approach. 
We note that the coupling strength 32 splits the com¬ 
mon energy E^ + fiwin = E2 + huj2 + fiioin of the states 
|3)at|?4, K)bos and |2)at|n, k-|- l)bos into an Autler-Townes 
doublet, and in such way a dressed picture is obtained, 
where the energies of the two dressed states are equal to 
Eh + huJin ± hg2i/K. If now we set hg2^/K = Eh — Ei 
we obtain for the lower energy of the dressed states 
Eh + huJin — hg2^/K = Ei + huJin. In such a way, for 
52^ giy/n, we obtain a one-dimensional coupled 
chain with a boson number gradient. This corresponds 
to the limit of degenerate qubit levels of the two-level 
quantum Rabi model, which was described in [ 23 , 1 ^ . 
This system is characterized with a strictly periodic be¬ 
havior of the populations and the boson number. We 
illustrate this on Fig.[ 3 l where we plot the time evolution 
of the populations of the three levels, the revival proba¬ 
bility |(4'(0)|4'(t))|^, and the boson number as a function 
of time. The initial condition is |4'(0)) = |l)at|0, «:)bos- 
The figure clearly reveals a periodic dynamics with a pe¬ 
riod of approximately 27 r/u;i. We also notice that the 
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FIG. 3. Time evolution of (a) the populations of states |l)at, 
|2)at and |3)at, (b) revival probability |('I>(0)||<I'(t))|^ and (c) 
boson number in the hrst mode. The energies of the bare 
states are Ei = 90&ui, E 2 = 0, E 3 = 100/kui. The coupling 
strengths are gi — l.Swi and % 2 %/k = E3 — Ei and the mode 
frequencies are uj 2 = lOOtui. 


mean boson number of the first mode never obtains too 
large values, which justifies the use of the approximation 
g 2 \/K ^ gi\/n. On the contrary, if the coupling strength 
is not tuned to E 3 — E 1 , we will get a non-periodic 
behavior of the dynamics. This is shown in Fig. (4] 


III. CONCLUSIONS AND DISCUSSION 

In this paper we have studied theoretically a mixed 
Rabi Jaynes-Cummings model of three-level atom inter¬ 
acting with two quantized bosonic modes, for which the 
RWA can be applied for one transitions solely. We have 
derived the differential equations, describing the evolu¬ 
tion of the system in Hilbert space, and shown that 
the weak-coupled bosonic mode can be exploited to con¬ 
trol the dynamics of the atom-field of the strong-coupled 
transition. In particular a transition from non-periodic 
to near-periodic behavior of the evolution has been found 


FIG. 4. Same as Fig. [3l but for hg 2 y/K = E 3 — Ei + 2huJi. 


by appropriate tuning of the control weak field. Such 
a transition has been related to the spectrum of the 
two-level quantum Rabi Hamiltonian, which shows an 
equally-spaced energy ladder structure in the limit of de¬ 
generate energy levels. As a final comment, let us briefly 
mention that the mixed Rabi Jaynes-Cummings model 
introduced in the present work can be simulated in an 
optical setting, following an approach similar to the one 
recently proposed and demonstrated in Refs. [H, [ 2 ^. 
The coupling scheme shown in Fig. [2] can be implemented 
in classical optics by using evanescently-coupled dielec¬ 
tric waveguides. For this aim, one should notice that 
Eqs. ([5]) are analogous to the coupled-mode equations de¬ 
scribing light transport in a semi-infinite one-dimensional 
photonic lattice in the tight-binding approximation with 
additional waveguides, which couple to only the even- 
number waveguides from the chain. The lattice shows 
a superimposed transverse index gradient and a non- 
uniform coupling constant between adjacent waveguides. 
The coupling constants are controlled by the distance 
dn between the waveguides and the index gradient wi is 
given by uji = 27rnsa/(i?A), where Ug is the refractive 
index of the substrate, A is the wavelength of light, R is 
the bending radius of curvature and a is the horizontal 
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spacing of the waveguides (see, for instance (311^). The rameters and can be experimentally determined, 
coupling strength J„ = hgiy/n between the neighboring 
waveguides is given to an excellent accuracy by the expo¬ 
nential law Jn = xexp(—adra), where x a are some ACKNOWLEDGMENTS 

constants which depend on the waveguide fabrication pa- 
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